Light-like Wilson Loops and Cusp Anomalous 
Dimensions in Nonconformal Gauge Theories 

Leopoldo A. Pando Zayas*, Daniel J. Phalen* and Cesar A. 

Terrero-Escalante^ 

66 '. 

f— ^ , * Michigan Center for Theoretical Physics 

f^*) ' Randall Laboratory of Physics, The University of Michigan 

^ ; Ann Arbor, MI 48109-1040 

''Centro de Investigacion y Estudios Avanzados 
Av. Politecnico Nacional 2508, Mexico D.F. 07360, Mexico 



o 



Abstract. We emphasize that nonconformal theories provide a natural playground for the ideas of 

the Maldacena conjecture opening the possibility of exploring properties that could potentially be in 

the same universality class as QCD. In particular, we discuss in detail how light-like Wilson loops, 

an important ingredient in the prescription for scattering amplitudes, can be described in a number 

of gravity duals of nonconformal gauge theories. We point out to a few universal properties and the 

Qh| prominent role of the strong scale. 
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The idea that the description of gauge theories might effectively involve strings has a 

long history starting from the dual models for hadronic resonances. In 1974 't Hooft 

Jll] showed that gauge theories admit a limit in which their perturbation theory can be 

00 ! thought of as coming from a string theory. Ten years ago, a precise formulation of this 

O ■ equivalence was put forward for a particular gauge theory. Namely, it was argued that 

>• ■ ^ = 4 supersymmetric Yang-Mills (SYM) is dual to type IIB string theory propagating 

k> ! on Ad S 5 X S^ ISSH]. In the ten years following the seminal work of Maldacena [2] the 

^ I gauge/gravity correspondence has seen a period of maturity. In its current incarnation 

the correspondence has turned into a very fruitful theoretical framework. 

Particularly exciting is the side of Maldacena's conjecture relating the strongly cou- 
pled sector of the gauge theory with weakly coupled gravity. Many problems in gravity 
have been solved during nearly a century of work. By turning strong coupling prob- 
lems of gauge theories into questions in classical gravity, the correspondence opened 
a window of opportunity into one of the long-standing puzzles of particle physics: the 
infrared properties of gauge theories. Some of the questions to be addressed here are 
confinement, the spectrum of hadrons and the properties of the strongly coupled quark- 
gluon plasma. However, important questions about nonprotected quantities lie beyond 
the domain of classical gravity and should properly be tackled using string theory. 

There are two main obstacles in applying the gauge/gravity correspondence to phe- 
nomenologically interesting theories. The high amount of supersymmetry and conformal 
symmetry is not present in nature. Although originally formulated as a duality between 



string theory on AdSs x S^ and maximally SYM in four dimensions, it has quickly be- 
come clear that such dualities extend to more general situations. Significant results in 
understanding QCD-like theories have been achieved by constructing supergravity the- 
ories whose dual gauge theories contain ^/T = 1, SU{N) supersymmetric Yang-Mills 
[|5i|6j. In a remarkable series of papers DTfl Klebanov and collaborators carried out a pro- 
gram that concluded with a supergravity background that is dual to ^ = 1 supersym- 
metric gauge theory that displays confinement and chiral symmetry breaking in the IR. 
This background is known as the warped deformed conifold or the Klebanov-Strassler 
(KS) solution. In the UV this background encodes a series of Seiberg duality cascades 
in the field theory. A fairly comprehensive description of how to encode ^ = 1 gauge 
dynamics into supergravity backgrounds was presented by Sonnenschein and Loewy 
in [|8|]. They analyzed the following properties: gauge group, supersymmetries, Wilson 
and 't Hooft loops and the corresponding quark anti-quark or monopole anti-monopole 
potentials, instantons and the U{1)r -^ Zin symmetry breaking, gluino condensation 
and spontaneous Z2N — > Z2 breaking, monopoles, domain walls, baryons and KK states. 
Moreover, a supergravity solution dual to the finite temperature warped conifold with 
cascading in the UV was found [9]. Evidence of a transition between this solution and 
the KS solution was reported in [10] (see also [11]). Such transition could be interpreted 
as the gravity dual to the confinement/deconfinement transition in QCD-like theories. 

Another important direction in which the Maldacena conjecture has been expanded 
is in going beyond the lowest energy states. Much of the developments in the con- 
formal version regarding beyond the supergravity approximation were extended to the 
nonconformal case. For example, the BMN limit \V2\ was also understood in noncon- 
formal situations; the analogous structure describes hadronic states for which an exact 
string Hamiltonian was presented [|l3ll . Using semiclassical quantization in the context 
of gauge/gravity correspondence, quantum corrections to Regge trajectories were calcu- 
lated [14], showing that quantum effects alter both the linearity of the trajectory and the 
vanishing classical intercept, J := a{t) = (Xq -\- a' t + /3 ^/F. 

More recently, Alday and Maldacena proposed a prescription to calculate the gluon 
scattering amplitudes for ^ = 4 using beyond the supergravity degrees of freedom 
IIS]. Further research in that direction has already gained momentum (for instances and 
references see the contribution by Alday and Maldacena to this volume). The logical 
next step is to apply their prescription to nonconformal situations. However, several 
obstacles arise. First of all, we have to deal with less symmetry, meaning that we need to 
use more powerful techniques to solve the corresponding problems. Fortunately, these 
problems reduce to systems of ordinary and partial differential equations, and we can 
use many analytical and numerical methods for studying and solving them. 

One goal of this contribution is to disp lay the power of these methods. Using as 
illustrative examples the Witten-QCD [fl6|] and KS L^] models, we are going to show 
that analysis of the light-like Wilson loop in the dual momentum space for nonconformal 
situations shows some universal properties. Based on the recognizable universality that 
we obtain, we expect that the analysis of gluon scattering in nonconformal theories 
would also yield valuable results potentially in the same universality class as QCD. 




FIGURE 1. Schematic representation of the Wilson loop. 



WILSON LOOPS IN GAUGE THEORIES 



The Wilson loop is a very important observable in gauge theories and in some ap- 
proaches it is treated as the central ingredient in the description of all gauge invariant 
quantities. Its description on the gravity side was proposed in [|l7l Iw] where it was 
suggested that the evaluation of the vacuum expectation value of a given Wilson loop 
on the contour ^ is given by the string worldsheet whose boundary coincides with ^: 
(WC^))=exp(-5(^)). 

The Wilson loop provided the first example of the correspondence that goes beyond 
the study of protected quantities at the supergravity level. Related developments allowed 
to understand that certain non-BPS operators can be described by classical string config- 
urations in a given supergravity background. For example, the large anomalous dimen- 
sion of operators of the form TrX^V/^^ . . . V^ ^X^, where X' are scalar fields in ^ = 4 
SYM, was computed in [19] by considering a large closed string spinning in AdS^, 
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where X = gy^N is the 't Hooft coupling. These operators are generalizations in ^ = 4 
SYM of the twist-two operators that play a central role in the deep inelastic scattering 
in QCD niol blll. They appear in the operator-product expansion of two electromagnetic 

currents Fa/xiV/i2 • • • '^a.-i'^"k. ^"^^ VYiii^^i ■ ■ ■ ^^.V- Soon after ^\, a description of 
these operators was found in terms of light-like Wilson loops in the AdS/CFT correspon- 
dence 1221 1231]. One important property of the anomalous dimensions of these operators 
is that it is conjectured to be in general of the form A — 5 = /(A) ln5, where /(A) is the 
cusp anomalous dimension and it is a function only of A. 

The Alday-Maldacena prescription for gluon scattering in ^ = 4 SYM |15|l involves 
computations of classical worldsheets in AdS^ x S^ that end on light-like segments. 
These configurations are light-like Wilson loops. It is crucial that the space where these 



Wilson loops live is the T-dual space to the original supergravity background. This con- 
struction motivates us to revisit and expand the analysis of worldsheet configurations 
ending on light-like segments in supergravity backgrounds dual to nonconformal theo- 
ries. Similar configurations have been considered before but with two main differences 
from those we present here. First, they were considered in the original supergravity back- 
ground, that is, without performing a T-duality and second the light-like segments were 
always placed in the asymptotic UV region. For example, the light-like Wilson loop on 



the Klebanov-Strassler background was analyzed in [|24|l . 



In this contribution we report on the situation where we perform a T-duality on the 

1 o 

R ' subspace where the field theory lives, and the light-like segments are placed in the 
region corresponding to the IR of the field theory. Since we work on the T-dual space and 
place the brane in an arbitrary position, the value of the cusp depends on the structure of 
the infrared region. This is a crucial difference with perturbative QCD where no access 
to the IR region is possible other than on dimensional grounds. 



Cusps anomalous dimension 

If a Wilson loop is evaluated over a closed contour that forms a light-like cusp with 



angle 7 in Minkowski space-time, then its expectation value is Il25ll . 



^,~y , (2) 

where A is a UV cutoff, m is an IR cutoff, and Tcusp measures the anomalous dimension 
of the Wilson loop. The cusp anomalous dimension is also an important observable in 
gauge theories. It controls the scaling behavior of various gauge invariant quantities like 
the logarithmic scaling of the anomalous dimension of higher-spin Wilson operators, 
double-log (Sudakov) asymptotics of elastic form factors in QCD, the gluon Regge 
trajectory, infrared asymptotics of on-shell scattering amplitudes, and it is important 
for resumming the effects of soft gluon emission in the study of QCD at colliders ll26ll . 
The cusp anomalous dimension depends only on the coupling constant and its expan- 
sion at weak coupling is known in QCD to three loops [27] and in ^ = 4 SYM to four 
loops [i28i1. I n the strong coupling regime of ^ = 4 SYM it was calculated to lowest 
order in rtl9l. l22Ll23|l . and to higher orders in [|29|,l30|]. In the strong coupling regime of 



cascading theories the first results were obtained in 112411 . For our results we needed to 
generalize this last approach by putting the brane in a different position. In the future 
we plan to calculate the scattering amplitudes and at that point we actually will require 
a surface instead of an infinitely extending cusp, that is, the surface for scattering ends 
on cusps and thus the surface does not extend to infinity. 

In this note we report mostly on work performed for WQCD and the KS backgrounds. 
However, for the strategy of understanding which precise properties of nonconformal 
theories are universal and which are model-dependent, it is important to discuss various 
cases. We are currently analyzing also the cases of the nonconfining cascading solution 
of ISD, the noncritical background R^'^ x5L(2,R)/[/(l) [31] and the noncritical solution 
with an infrared fixed point [32]. 



The general set up for the cusp 

Let us consider a general metric of the form: ds^ = q^{r)dx^dxpi + p^{r)dr^. We 
perform a T-duality along the x^ coordinates and arrive at a metric of the form 

ds^ = -^dy^'dy^+p\r)dr\ (3) 

We now consider embedding a string world sheet into this metric. For world sheet 
coordinates (t, a), this embedding is yo = e^cosha, yi = e^sinha, r = r(T, a) . The 
action is thus 

S = ^,jdadx'^^\-pWe-^'{{d,rY-{d,rY). (4) 

The methodology is then to derive the equation of motion for r, to solve it, substitute 
the solution into (lH) and integrate for the action. This integral typically diverges, so an 
UV cutoff must be chosen as upper limit of integration. Given that S = \nW , the cusp 
anomalous dimension can be obtained from ^ as. 

The cusp in^ = 4 SYM 

Here we present elements of a systematic analysis to the derivation of the solution 
for this most symmetric case. It might seem like an overkill, however, reproducing 
analytical results for simple cases helps to adjust the numerical methods for handling 
less symmetric situations. Here we study a building block for the configuration that 
leads to the scattering amplitude. 

Let us consider the cusp as described in nlSn . For AdS^ we have that in dH) ^^ = r^/R^ 
and p^ = R^/r^. Moreover, if we consider an ansatz of the type r = e^w(T), then we find 
that the action is: 

Following the methodology listed at the end of the previous subsection we derive the 
corresponding equation of motion which, for the purpose of numerical integration, is 
better to write as the following system of first order ordinary differential equations, 

vv = V, V = — (2 — 4vw — 2v — 3w +w +VW — V w — V w) . (7) 

w ' 

Here dot stands for derivative with respect to T. We solved this system using the seventh- 
eight order continuous Runge-Kutta method. Thanks to its adaptive scheme, this method 
provides a great control upon the output accuracy. Typical solutions are plotted in figl2l 
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FIGURE 2. Left: The flow defined by system ^. Center: The flow near the stable node. Right: The 
flow near the saddle point. 



An important point here is that numerical solutions are not necessarily black boxes. 
Qualitative and asymptotic analysis help to understand the numbers which in turn might 
guide the analysis. For instance, in this case we can use the theory of dynamical systems 
to understand the contents of figures [2l This theory includes a set of general techniques 
for the analysis of phase spaces that allows us to understand better the asymptotics of 
solutions and their dependence on the initial conditions. This is a natural way of tackling 
problems about classical solutions in string theory; we are not aware of previous works 
emphasizing these systematic methods; most of the literature emphasizes exploiting 
symmetries. While of paramount importance in AdSs x S^, most symmetries are lacking 
in the context of supergravity backgrounds dual to nonconformal theories. 

Setting the vector field to zero in system (|7]), we see that there are two fixed points. 
They are respectively located at (1,0) and (-\/2, 0). If these singular points were hyper- 
bolic, then the Hartman-Grobman theorem [33] ensures that, in a small neighborhood of 
each hyperbolic fixed point, the flow defined by the full nonlinear system (|7]) is topologi- 
cally equivalent to its linearized version. A singular point is hyperbolic if the eigenvalues 
of the Jacobian of the linearization around the point have non-zero real parts. For the first 
point we obtained Ai = — 1 and Ai = —2. This indicates that in its neighbourhood the 
flow behaves like near a radial sink (stable node), as it is shown in the center of figEl 
For the second singular point Ai = — 1 -f ^/3 and Ai = — 1 — ^/3. So, in its neighborhood 
the flow behaves like near a saddle point, as shown to the right in figl2l 

Having tested that the numerical and analytical results coincide, we can move forward 
and integrate Q numerically. Taking A = r,nax, we evaluate the action with a set of 
different upper limits t,nax which are defined by rmax- Then using ^ = j^-^f—, we 
estimate the cusp anomalous dimension using ([5]). The outcome for this case nicely 
replicates the result in [15]. Because of lack of space, we refrain of presenting it here, 
and move forward to the nonconformal cases. 



CUSP ANOMALOUS DIMENSIONS IN NONCONFORMAL 

GAUGE THEORIES 

Witten QCD 

In this subsection we calculate the cusp anomalous dimension for the Witten back- 
ground [16]. With that aim we rewrite the metric in the relevant IR regime. The ten- 
dimensional string frame metric and dilaton of this model are 

K yUQ U T\U) 

1/3 , 1 ,,3/4 

/(") = 1-^' R={KNgs)^a'\ ^'^ = Ss-^,. (8) 

The geometry consists of a warped, flat 4-d part, a radial direction u, a circle parame- 
terized by 9 with radius vanishing at the horizon u = uq, and a four-sphere whose volume 
is instead everywhere non-zero. It is non-singular at w = uq. In the m ^ oo limit the dila- 
ton diverges: this implies that in this limit the completion of the present IIA model has 
to be found in M-theory. The background is completed by a constant four-form field 
strength F4 = 37?-^ 0)4, where (04 is the volume form of the transverse S^. 

The main gauge theory parameter we will use in the following is the KK mass scale 
l/i?0 — 3mo/2, where 11^ = uq/R? . As can be read from the metric, mo is also the typical 
glueball mass scale, and its square is proportional to the ratio between the confining 
string tension Tqcd and the UV 't Hooft coupling A. The supergravity approximation 
works in the regime opposite to that in which the KK degrees of freedom decouple from 
the low energy dynamics. Condition Tqcd ^ "^0 implies A ^ 1, which is beyond the 
supergravity regime of validity. 

Assuming that u{x,o) = r(T), with q{r)^ = (u/R)^^^, p{r)^ = {R/uf^^ x 1//, and 
introducing the new coordinate t = e^, so that e^'^d /dx = d/dt, the action ^ is now: 

R^/^ r r 1 / 

'da dtt^j^J\-{d,rY/f. (9) 



2na' J J r>l^ 



A convenient way to write the equation of motion is: 

1 [2rV3-3/'r3(r3-2r3)r2-2r4(r3-r3)2r + 3/'(r3-r3)2] 



2 ^4^(r3_^3) 



(10) 



where r and t are now measured in units of R, and a dot stands for derivative with 
respect to t. Note that r = const = ro is not properly a fixed point of the action. A typical 
solution of this second order non-linear non-autonomous equation is presented in the 
left of figS As with the rest of the examples plotted in this subsection, it was obtained 
by setting ro = 1, ?o = 10"^, '"(0) = 1 + 10"^ and r(0) = 0. To understand this result let 
us analyze the asymptotic regimes. 
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FIGURE 3. Left: a typical solution for r in the WQCD model. Center: error of the approximate solution 
(fT2] i compared to the numerical solution. Right: phase space for very large t. Dr stands for the derivative 
of r with respect to t. 



Asymptotic behavior of the solution 



First, we look for a solution for {t 
(fTOl) the relevant terms are, 



-to) ^0 and e = r(0) — ro^O. Balancing equation 



3 -3 



-ir'-rl) 



t r^ 



(11) 



For this balance to be consistent it is required the absolute value of the second derivative 
to be greater than the neglected terms. It can be readily tested that this condition is 
fulfilled in the limits of £ and t we are considering. 

We inmediately note that any solution to (fTTI) with initial condition r = at ?o 7^ is 
a constant solution. A solution to this equation given as an expansion in Taylor series is. 



r={ro + £)+£{t-to) + 



e{r^£ + Srge^ + 3rQ£^ + £'^- Srg - 3ro£ - e^ 
2to{3rl + 3ro£ + £^) 



(12) 



It can be seen that, in the limit e — > 0, ro actually acts as a fixed point of the system. We 
tested this result numerically. In the center of figl3]we can see that the difference between 
(fT2)) and a corresponding numerical solution is negligible. Thus, the approximation 
works in the range of t and e considered. 

Let us now consider the asymptotic behavior of the solution for t ^ ^ Now the 
equation of motion (flOl) becomes 



'+2 



3r^ir^-2rl) 



(r3 



fo)r 



3r3 

2~ 



0. 



(13) 



With regards to the consistency of this balance, from figl3]it can be assumed that the large 
t behavior of the solutions is linear. Then, the second derivative goes to zero and that is 
also true for the neglected terms, {Pr^)/{t{r^ — Tq)) — r/t . For further verification we 
again compare the analytical and numerical results. Note now that the second derivative 
is zero when the first derivative is ±{P — rQ){{r^ — 2rQr)^'^r) whose limit as r — > oo 





FIGURE 4. Left: fit to 1 jt of the argument of the squared root in the WQCD action. Center: plot of the 
integral as function Qitmax- Right: fit of C as function of z. 



gives ±1. So, for any ro, some solutions converge to a linear solution with slope 1 and 
some diverge from a linear solution with slope —1. In the right of figl3]the phase-space 
of (fT3l) is plotted. We see that, indeed, for any ro, r = 1 is an attractive set and that the 
solutions diverge from r = — 1 . Now, recall that we are analyzing the behavior for r ^ oo 
and, according to the analysis of the solutions when Z^ ^ 0, the interesting us region is 
r > ro and r > 0. It means that the interesting us asymptotic behavior is given by the 
linear solution with slope equal to 1 . 

Summarizing, the analysis and the numerical results for the range of ro tested (ro G 
[10~^, 10] and fo = 10^^) indicate that the generic behavior of the solutions is as follows: 
ro acts as a fixed point of this system; with an initial condition near ro, the solution starts 
with a regime where r ?^ ro. Next there is a transient regime where the solution smoothly 
changes until it finally sets into a linear behavior with unit slope. The closer the initial 
condition to the horizon, the larger the length of the initial regime and the shorter the 
transition regime. The length of the time intervals are to be considered relatively, since 
the figures show always the same behavior if the right lengths of the t interval are chosen. 



The cusp anomalous dimension 



Next, following our methodology, we evaluate the numerical solution in the action 
and integrate it numerically. To accurately compute the actions we used an adaptive 3-5 
Simpson's quadrature. If numerical factors are not taken into account, the integrand is, 

t I r^ii — f^^—l 

-jj2 y .3_. — . We have already seen that r ^>- 1 while r ^ oo. Therefore, if we expect 
this integral to diverge logarithmically, we need the argument of the squared root to 
behave like l/t. Indeed, using a best fit to the numerical results, we have verified that 
this is the case, and we present an example in the left side of fig HI In the center a typical 



result for the action as a function of t„ 



°^ r„ 



A is plotted. 



The coefficient of dS/dlnr,nax is A^/C, where A stands for the numerical factors in 
front of the integral, C = ro{bz^ + ^o + biz + b2Z^ + hz'^), z = r{0)/ro - 1 and p = 
0.27895, Z7 = 0.436898, Z70 = 0.016228, Z^i = 0.359128,^72 = -0.026669,^73 = 0.01103. 
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FIGURE 5. There is a fixed point of the Klebanov-Strassler equations of motion at po = 0. This plot is 
for e = 10^1 



The numbers were obtained by the best fit C{z) shown in the right of fig® 



The warped deformed conifold 

The Klebanov-Strassler background ^^ is very well known in the literature of gravity 
duals to confining gauge theories, a full presentation of its metrics and action, suited for 
numerical calculations can be found in the appendix A of [l(j]. 

Using coordinates where t = e'^ and setting p{t, o) = p{t), the equation of motion for 
this background is 



£^l^h'l^{p)p 



6tK2{p)Jl 



£4/3 . 2 



(14) 



th^/^ip). 



£4/3 



( dpKp) 



e^/^p%Kip) 



et^K^ipf \^2h{p) t^K\p)il 



£4/3 




This equation was solved by setting p (0) = £ ^ and p'(0) = 0, and then using a fourth 
order Runge-Kutta method for the evolution of the solution as a function oft = e^. The 
numerical integration necessary for evaluating the warp factor h{p) was done using a 
composite Simpson method. A typical solution is presented in figlS] An analysis similar 
to the previous subsections shows that there is a fixed point at p = 0. The solution 
initially stays close to the fixed point, then there is a transition to a logarithmic looking 
function. The region close to the fixed point is larger depending on how close to the fixed 
point the solution begins. 

To show that the solution grows logarithmically, we examine the form of the metric for 
large t. In this regime, p is also getting large, so we have K{p) — > B exp(— p/3), where 



5 is a constant. The KS metric now has the form ds^ = h ^'^{p)dx^ + h^'^e^P'^dp^. 
Substituting r = e^'^ and peri'orming a T-duality on the x coordinates, the metric is 

ds^ = h^/^{3logr){dx^ + dr^). (15) 

Now, in the large r (large p) approximation h{p) -^ A(31ogp — l/4)/r^, where A is a 
constant, and one gets 

J/ = ^3logp-l/4 '^''^Y'^ (16) 

which is AdSs with a logarithmic prefactor. To a first approximation, we can ignore the 
prefactor which will change slowly, to the metric, to see that the cusp solution of AdS^ 
which is r = constant x t will hold, which means that p ~ log? for large t. 
This way, the result for KS resembles very much the one obtained for WQCD. 



DISCUSSION AND OUTLOOK 

In this contribution we have shown that light-like Wilson loops in nonconformal theories 
display some rather universal features. We hope that this is the kind of question that can 
be asked and whose answer should be trusted to be in the same universality class as 
QCD. The light-like Wilson loops we have considered here are building blocks of the 
scattering amplitudes. We hope, as is the case in the perturbative treatment of QCD 
amplitudes, that the main contribution to scattering amplitudes will come from the cusp 
region and therefore, our calculations will be a crucial building block in the amplitudes. 
It is important to note that, although not included in this contribution, we have 
performed similar analysis for various theories and the behavior is rather similar to 
the one observed in the two examples considered explicitly here. Basically, there is a 
natural infrared regulator which is characterized by the distance between the position 
of the brane and the characteristic strong scale which in most supergravity solutions is 
given by the region where a cycle collapses (ro for WQCD and p = for KS). We want 
to emphasize that in the perturbative approach to gauge theories the infrared regulator 
normally used is formal and has nothing to do directly with the strong scale. We thus 
hope that our strong coupling calculation will help shed some light into the structure of 
IR singularities. 
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